WMAP, neutrino degeneracy and non-Gaussianity constraints on isocurvature 
perturbations in the curvaton model of inflation 
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In the curvaton model of inflation, where a second scalar field, the 'curvaton', is responsible for 
the observed inhomogeneity, a non-zero neutrino degeneracy may lead to a characteristic pattern of 
isocurvature perturbations in the neutrino, cold dark matter and baryon components. We find the 
current data can only place upper limits on the level of isocurvature perturbations. These can be 
translated into upper limits on the neutrino degeneracy parameter. In the case that lepton number 
is created before curvaton decay, we find that the limit on the neutrino degeneracy parameter is 
comparable with that obtained from Big-bang nucleosynthesis. For the case that lepton number is 
created by curvaton decay we find that the absolute value of the non-Gaussianity parameter, |/ n i|, 
must be less than 10 (95% confidence interval). 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

In single field models of inflation, the inflaton gives rise 
to an era of accelerated expansion in the early universe 
during which quantum fluctuations in the inflaton field 
are expanded beyond the horizon size [1]. At the end 
of inflation, the inflaton decays into the standard matter 
fields. The primordial fluctuations, imprinted onto the 
gravitational potential in form of small perturbations, be- 
come the "seeds" that start the formation of large scale 
structure. As there is only a single degree of freedom, 
the inflaton's value, the resulting perturbations are al- 
ways adiabatic, i.e. the ratio of the number densities of 
different particle species is spatially homogeneous, there 
is no isocurvature mode. 

Recently a different model has been proposed, the cur- 
vaton scenario [2, 3], which has already received consider- 
able attention [4-7] . In this model the accelerated expan- 
sion is still provided by the inflaton, but the primordial 
perturbations are generated by another scalar field, the 
curvaton. The curvaton energy density is small compared 
to the inflaton energy density and being a light field, it 
acquires an almost scale-invariant spectrum of perturba- 
tions during inflation. After the decay of the inflaton, the 
isocurvature perturbation in the curvaton fluid is trans- 
formed into an adiabatic perturbation [5, 8, 9] and the 
curvaton decays into radiation and other particles. A 
similar model to the curvaton scenario, in which a sec- 
ond scalar field determines the inflaton decay rate, has 
also been proposed recently [10]. 

The curvaton scenario requires different constraints on 
the inflationary potential to standard inflation and might 
therefore be better suited for inflation at a lower energy 
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scale [11, 12]. It also allows for the possibility of an 
isocurvature mode in the standard matter fields [2, 3], i.e. 
the ratio of particle number densities of different parti- 
cle species is not necessarily spatially homogeneous. This 
arises since the final constituents of the standard big bang 
model (photons, neutrinos, baryons and cold dark mat- 
ter (CDM)) do not all necessarily originate entirely from 
the spatially inhomogeneous curvaton, but may originate 
instead from the spatially homogeneous inflaton or some 
other source. 

In Ref. [7] the isocurvature components in the CDM 
and baryons that could arise in the curvaton scenario 
were constrained using WMAP and other observational 
data. It was found that the current data rule out either 
CDM or baryons being created before curvaton decay. In 
this article we examine what the current observational 
constraints are on neutrino isocurvature perturbations 
created by curvaton decay. These can only arise if there 
is a non-zero neutrino degeneracy parameter £ = n v /T v , 
where [i v is the neutrino chemical potential and T v is the 
neutrino temperature [3]. There may be a separate de- 
generacy parameter for each of the three species of neutri- 
nos, but the recent atmospheric and solar neutrino data 
indicate that approximate flavor equilibrium between all 
active neutrino species is established well before the BBN 
epoch [13] and hence we take the degeneracy parameters 
to be equal. 

Various studies have constrained the value of £ using 
CMB and other large scale structure data in the adia- 
batic case [14]. However, the tightest constraints still 
come from Big-bang nucleosynthesis (BBN) [15]. In this 
article we examine how the constraints on £ change in the 
curvaton scenario where neutrino isocurvature perturba- 
tions can be induced by a non-zero £. It turns out that 
the neutrino isocurvature perturbations then also induce 
baryon and CDM isocurvature perturbations [3]. It is 
therefore not possible to use constraints that have simply 
considered the addition of a neutrino isocurvature mode 
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such as Ref. [16] or those that have considered arbitrary 
combinations of different isocurvature perturbations [17]. 

In the next section we evaluate the possible magni- 
tudes of the neutrino isocurvature perturbations. Then 
in Sec. Ill we estimate analytically the effects of the 
isocurvature perturbations on the temperature fluctua- 
tions of the CMB and their dependence on £. A full 
likelihood analysis of the allowed ranges of £ is given in 
Sec. IV. Finally, the conclusions and the discussion are 
given in Sec. V. 



II. CURVATON GENERATED NEUTRINO 
ISOCURVATURE PERTURBATIONS 

A useful gauge-invariant variable for characterizing in- 
homogeneity on large scales is [18-20] 



P 



(1) 



where p and H are respectively the total energy density 
and the Hubble parameter, the dot denotes differentia- 
tion with respect to coordinate time, and the perturba- 
tion in a background quantity x is denoted by Sx. For 
a critical density Universe, tp is related to the intrinsic 
curvature of a spatial hypersurface by [20] 



(3) i?= ± V 2 V)) 



(2) 



where a is the scale factor and Vj is the covariant deriva- 
tive. It is possible to choose the coordinate system such 
that certain combinations of perturbation variables are 
set to zero. This is known as choosing the gauge or frame 
[20-22] . In this article we will work in the flat gauge [20- 
22] in which ip — 0, and Eq. (1) becomes 
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where from now on Sx refers to the perturbation in x 
when the flat gauge is chosen. The energy conservation 
equation is given by 



p + 3H (p + P) = , 



(4) 



where P is the total pressure. After curvaton (a) decay, 
the component fluids are photons (7) , neutrinos (u) , cold 
dark matter (cdm) and baryons (B). The total energy 
density and pressure are related to density and pressure 
of a component i by p = pi and P — Pi . The com- 
ponent pressures are given in terms of the corresponding 
component densities: Pi = pi/ 3 for photons and neutri- 
nos and Pi = for the curvaton (after the end of infla- 
tion), CDM and baryons. 

We assume a critical density Universe and that the 
neutrinos are effectively massless. 

The isocurvature perturbations are given by [5] 



where i is one of a, v, cdm or B and 

Pt 



(6) 



If there is no energy transfer to and from a fluid so that 
it satisfies the energy conservation equation: 

p\ + 3H( Pl + Pi) = 0, (7) 

then its Q is conserved on large scales, Q = [20]. 

The occupation number for a neutrino species (i = e, 
p or t) with energy E is given by 

f i (E) = {exp(E/T v Tti) + l}- 1 , (8) 

where T v is the neutrino temperature and the minus is 
for neutrinos and the plus for anti-neutrinos. The degen- 
eracy parameter is & = pijT v where pi is the chemical 
potential of species i. However, using the Large Mixing 
Angle (LMA) solution leads to the degeneracy parameter 
being approximately the same for all three species [13], 
i.e. £j = £. After positron annihilation, the asymmetry 
parameter, £, is constant in time [3]. When £ is non-zero, 
the difference between the number densities of neutrinos 
and anti-neutrinos is non-zero. We denote this quantity 
riL as this difference is equal to the lepton number. Given 
the constraints on charge asymmetry, any non-negligible 
lepton number must be due to the neutrinos. There can 
be no neutrino isocurvature perturbation when til = 
as then the neutrino energy density is determined solely 
by the photon energy density before neutrino decoupling 
[3]. 

Analogously to the case of energy density, Eq. (6), the 
inhomogeneity of lepton number can be characterized by 



(9) 



where as specified before we are using Shl in the flat 
gauge. If the lepton number density is conserved 



n L + 3Hn L = , 



(10) 



then (l = on large scales [3, 6]. 

The isocurvature perturbation of the lepton number 
perturbation is then given by 



5 L = 3(Cl - C 7 ) , 



(11) 



which is related to the neutrino isocurvature perturbation 

by [3] 



where 



3(C* - Cy) : 



(5) 



B = 



45 B 2 
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(13) 
(14) 
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and 5' = 1+3(£/tt) 2 . The 3.04/3 term in the definition of 
A takes into account the non-equilibrium heating of the 
neutrino fluid and finite temperature QED corrections 
[15, 23, 24]. 

If the lepton number is created well before the curvaton 
decay then it will have negligible perturbations (Cl = 0) 
and so from Eq. (11) 



-3C 7 . 



(15) 



We will not make the assumption |£| -C 1, which was used 
in Ref. [2] as we will be also considering the case where 
£ may be large. To evaluate the effect of S u on obser- 
vations we need to express it in terms of the adiabatic 
perturbation (. To do this we can use Eqs. (3), (4) and 
(6) and that only the photons and neutrinos contribute 
significantly to the density in the primordial era to get 



C = (1 - Ru)C-y + RuCu ■ 



where 



R v — 



Pv + Py 

which can be evaluated using [15, 23] 



Pv 
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Combining Eq. (16) with Eq. (5) for i — v we get 
C = C7 + ~zRvSv 1 



(16) 



(17) 



(18) 



(19) 



then solving Eq. (19) for £ 7 and substituting into Eq. (15) 
gives 



-3( + R V S V . 



(20) 



Substituting Eq. (20) into Eq. (12) and solving for the 
neutrino isocurvature perturbation gives 



S v = 



135B 2 



A5B 2 R V - IB' A 



(21) 



If the lepton number is created by curvaton decay then 
[3] 



Cl = C. 



c 



(22) 



where using the sudden decay approximation, r w 
Pa/{p<y + Pi) evaluated at the time of curvaton decay. 
The approximation is accurate up to about 10 % error 
even when the decay is not sudden [5]. Using Eq. (5) 
(with i = v), Eqs. (11), (16) and (22) gives 



1 - r 

Sl = 3 C + RvS v ■ 



(23) 



Substituting Eq. (23) into Eq. (12) and solving for S v we 
get 



S v = 



r - 1 



135B 2 



r A5B 2 R U - 7B'A 



1 A £ ' 



(24) 



which is valid even when £ 2 is large. As can be seen, the 
neutrino isocurvature perturbation resulting from lepton 
number created by curvaton decay, Eq. (24), differs from 
the neutrino isocurvature perturbation result when lep- 
ton number is created before curvaton decay, Eq. (21), 
by a factor which is only a function of r. 

The final possibility is for the lepton number to be 
created after curvaton decay which corresponds to r = 
1 in Eq. (24) and so in this case there is no neutrino 
isocurvature perturbation. 

If there is a neutrino isocurvature perturbation it in- 
duces a CDM and baryon isocurvature perturbation [3]. 
In the case where CDM or baryons are created after cur- 
vaton decay [3] : 



<5 c dm/B — Ri>S v , 



(25) 



where 5 c( j m /B stands for S ci \ra or 5b as the effect is the 
same on either one. When CDM/B is created before cur- 
vaton decay there is an additional residual contribution 
[3]: 

Scdm/B = RuS u - 3C • (26) 

When CDM/B is created by curvaton decay 

1 - r 

Scdm/B—RvSv+Z (. (27) 

r 

However, to concentrate on the role of neutrino isocurva- 
ture perturbations we will only consider the case where 
baryon number and CDM are created after curvaton de- 
cay, Eq. (25). 



III. NEUTRINO DEGENERACY 
MODIFICATIONS OF THE SACHS- WOLFE 
EFFECT 

In this section we examine analytically the effect of a 
neutrino degeneracy in the curvaton scenario. 

In the presence of isocurvature perturbations the tem- 
perature fluctuations due to the Sachs- Wolfe effect [25, 
26] on large scales are given by [7] 



AT 



C' a ~ 7 {RbSr + J?cdm«5cdm) + —RvS v , 



where C rad is the value of £ in the radiation era and 

PB 



Rb — 

Rcdm 



PB + Pcdn 
Pcdm 

PB + Pcdu 



(28) 

(29) 
(30) 



Substituting the induced baryon isocurvature perturba- 
tions Eq. (25) into Eq. (28) and making use of Eqs. (29) 
and (30) gives 



(31) 
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The current 95% confidence interval (CI) limits on £ from 
BBN arc [151 



-0.03 < £ < 0.11. 



(32) 



So in order to determine whether within these limits the 
isocurvature perturbations can still play a role we can use 
£ <C 1 which gives from Eqs. (14), (17) and (18) gives 



R u ps 0.41. 
Substituting Eq. (33) into Eq. (31) gives 
AT 1 / irad 1 



T 



C raa + 



(33) 



(34) 



which shows that the neutrino isocurvature perturbation 
has a similar sized effect as the adiabatic perturbation on 
large scales. 

For the lepton number created before curvaton decay 
with £ < 1, Eq. (21) gives 

which when substituted into Eq. (34) results in 

^«i(-l + £ 2 )C rad . (36) 

Hence in this case the neutrino degeneracy decreases the 
temperature fluctuations on large scales. However, for 
the BBN limits of |£| < 0.11 it follows from Eq. (36) that 
the isocurvature perturbation has a negligible effect on 
the temperature fluctuations when the lepton number is 
created before curvaton decay. 

For lepton number created by curvaton decay with £ -C 
1, Eq. (24) gives 



l-rl35 
~r T 



^From Eq. (37) and Eq. (34) we then find 
AT _ 1 ( r ~ 1 



-1 



-en 



rad 



(37) 



(38) 



As can be seen from Eq. (38), in the case where the lepton 
number is created by curvaton decay, a non-zero £ adds 
to the magnitude of the fluctuations on large scales as 
< r < 1. Limits on r can be placed using the limits on 
the non-Gaussianity parameter f n \ [3, 27] 
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(39) 



The current WMAP limits two sigma confidence interval 
is [28] 



-58 < /„i < 134 . 



(40) 



Inverting Eq. (39) this leads to r > 0.02 at two sigma 
confidence level. So from Eq. (38), it appears that there 
may be scope for a large contribution by the isocurvature 
term. We will test this with a full likelihood approach in 
the next Section. 
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FIG. 1: An illustration of the effect of the lepton degeneracy 
(£) on the spectrum of fluctuations. The cosmological param- 
eters (including the amplitude) for all the spectra are taken 
from the maximum likelihood when £ = 0. The solid line 
is a plot of the pure adiabatic spectrum. The short dashed 
fine is for the case when the lepton number is created before 
curvaton decay and £ = 0.3. The long dashed line is for the 
case when the lepton number is created by curvaton decay 
with £ = 0.3 and the fraction of the density made up by the 
curvaton at curvaton decay (r) is set to 0.5. The dotted line 
is the spectrum for when the lepton number is created after 
curvaton decay (so that there are no isocurvature modes) and 
£ = 1. The binned WMAP TT data is also plotted. 



IV. DATA ANALYSIS 

The effect of a non-zero £ on all cosmological scales can 
be obtained by inputting the isocurvature perturbations 
derived in Eqs. (21) and (24) and the induced baryon and 
CDM isocurvature perturbations, Eq. (25), into CAMB 1 . 
To account for the effect on the background dynamics the 
number of massless neutrinos needs to be set to 



N v = 3A 



(41) 



where A is defined in Eq. (14). Some example cases are 
plotted in Fig. 1. This figure confirms the direction of the 
effect evaluated analytically on large scales in Eqs. (36) 
and (38), i.e. a non-zero £ decreases the power on large 
scales when the lepton number is created before curvaton 
decay and increases it when the lepton number is created 
by curvaton decay. As can be see from Fig. 1, there is 
a similar effect on the cxtrcma of the acoustic peaks but 
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not on the rise and fall of the acoustic peaks. There is 
also a slight shift in the phase of the oscillations. As can 
also be seen in Fig. 1, a similar effect of a non-zero £ 
occurs when there are no isocurvature modes [14, 29] but 
then there is no effect on the Sachs- Wolfe plateau. Note 
that the decrease or increase is fairly constant over large 
scales and so is not suitable for resolving the quadrupole 
problem [30]. 

The current BBN constraints on £ are given in Eq. (32) . 
They rely on using the large angle solution in which all 
the £ for the three neutrino species become equal [13]. 
The actual constraint is on £ e , the electron neutrino de- 
generacy, due to its effect on the neutron to proton ratio 
prior to BBN. There have been attempts to check the 
BBN constraint by independently using CMB and other 
large scale structure data [14]. This has been done in the 
context of adiabatic perturbations and so the effect is due 
to the increase of density caused by an greater effective 
number of neutrinos, Eq. (41). 

In this Section we test what the constraints are when 
lepton number is created before and by curvaton decay. 
The adiabatic case is equivalent to when lepton number is 
created after curvaton decay as it corresponds to setting 
r = 1 in Eq. (24). A modified version of COSMOMC 2 
[31] was used for the likelihood analysis which included 
WMAP temperature and temperature-polarization cross- 
correlation anisotropy [32], ACBAR [33], CBI [34], 2dF 
galaxy redshift survey [35] and Hubble Space Telescope 
(HST) Key Project [36] data. The Universe was taken 
as flat and the neutrinos as massless. Reionization was 
parameterized using optical depth [37]. Each simulation 
used five separate Markov chains and the convergence 
criteria was taken to be the variance of the chain mean 
divided by the mean of the chain variance to be less than 
0.2 for each parameter. 

As COSMOMC produces Monte Carlo samples of the 
output variables, marginalization simply entails extract- 
ing the samples of the variables of interest from the multi- 
dimensional samples. These samples can then be used to 
evaluate confidence intervals or to give estimates of the 
probability distribution by using histograms. Functions 
of the variables can also be analyzed just by taking func- 
tions of the corresponding samples. 

The different scenarios that were checked were: lepton 
number created before curvaton decay, Eq. (21), lepton 
number created by curvaton decay, Eq. (24), and lepton 
number created after curvaton decay. COSMOMC was 
modified to take extra variables £ and r and their effect 
on the isocurvature modes, Eqs. (21), (24) and (25), and 
effective number of neutrinos, Eq. (41), were also added. 
Modifications were also made to allow the inclusion of 
a BBN prior on £ which is a Gaussian distribution with 
mean 0.04 and standard deviation 0.035 [15]. The prior 
for / n i was a Gaussian distribution with a mean of 38 
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5,/C rad 
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BBN only 


0.1 






before curvaton decay 


0.2 


-0.08 




BBN + before curvaton decay 


0.1 


-0.02 




by curvaton decay 


0.8 


0.07 


-4 


BBN + by curvaton decay 


0.1 


0.04 


-10 


after curvaton decay 


1.6 
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TABLE I: Upper bound (95% CI) constraints on the mag- 
nitudes of the lepton degeneracy (£), the neutrino isocurva- 
ture perturbation (<S„) and the non-Gaussianity parameter 
/ni- The 'BBN' label means the constraint from BBN using 
the large mixing angle solution [15] is included. The con- 
straints from WMAP, ACBAR, CBI, 2dF and HST on the 
various curvaton mechanisms of generating the lepton num- 
ber are given. 



and standard deviation of 48 [28]. This then imposed a 
prior on r through Eq. (39). 

In the cases where the BBN constraint was not used, 
only positive values of £ were sampled as then the results 
were even in £ as Eqs. (21), (24) and (41) are even in £. 

The constraints on £ for the various possibilities consid- 
ered are given in Table I. The marginalized distributions 
for £ are also plotted in Fig. 2. 

The BBN constraint is still tighter than that obtained 
by using the CMB, 2dF and HST data on any of the 
curvaton scenarios. 

The case where lepton number is created before curva- 
ton decay is only twice as broad without the BBN con- 
straint. Also, in this case, the isocurvature contribution 
is almost 10% as large as the adiabatic contribution but 
this gets dramatically reduced when the BBN constraint 
is added. 

The case when the lepton number is created by curva- 
ton decay is less constraining on £ even though the data 
do allow about as much positive isocurvature perturba- 
tions as negative. This is because there is an additional 
degeneracy with r as seen in Eq. (24) and so also in the 
two dimensional marginalized distribution of £ and r as 
seen in Fig. 3. Adding the BBN constraint also dramati- 
cally decreases the amount of isocurvature perturbations 
in this case as the smaller values of r needed to com- 
pensate for the lower £ is disfavored by the WMAP non- 
Gaussianity constraints on f n \. However, smaller values 
of r are more likely with the BBN constraint and these 
lead to values of f n \ large enough in magnitude to pos- 
sibly be detectable by the Planck satellite which may be 
sensitive to |/ n i| > 5 [38]. 

In the case where lepton number is created after cur- 
vaton decay, the constraints on £ are the broadest. In 
this case there are no isocurvature perturbations and a 
non-zero £ only effects the background dynamics through 
modifying the effective neutrino number, Eq. (41). This 
case has been studied before with similar results [14]. 
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-2-10 1 2 

FIG. 2: Comparison of marginalized probability distributions, 
using CMB, 2dF and HST data, for lepton degeneracy (£) 
when lepton number is created by curvaton decay (short- 
dashes), for lepton number created after curvaton decay (dot- 
ted) and for lepton number created before curvaton decay 
(long dashes). The BBN constraint on £ is also plotted (solid 
line). 



V. CONCLUSIONS 

We have examined the current cosmological constraints 
on the neutrino degeneracy £, in the curvaton scenario, 
using WMAP, ACBAR, CBI, 2dF, HST and BBN data. 
The expressions for the isocurvature perturbations were 
evaluated to all orders in £ and the effects of the in- 
duced baryon and CDM isocurvature perturbations were 
accounted for. Analytical formulas were given for the 
modification of the Sachs- Wolfe effect on large scales. 

Comparisons were made with the BBN constraint 
which results from the effect of the neutrino degeneracy 
on the pre-BBN proton to neutron ratio and using the 
large mixing angle solution in which the neutrino degen- 
eracies of the different species take on the same value. 
It was found that the constraint on £ was comparable 
to that from BBN when the lepton number was created 
before curvaton decay. As seen in Fig. 1, a non-zero £ 
affects the acoustic peaks whose data will dramatically 
improve in the future. So the CMB constraints on £ 
may also improve dramatically. Without the BBN con- 
straint, the neutrino isocurvature perturbation was found 
to be smaller than 10% the size of the adiabatic mode 
(95% CI), while including the BBN constraint reduced 
this by a factor of a few. In the case where the lepton 
number is created by curvaton decay, there is a degen- 
eracy between the value of the curvaton energy density 
at the time of curvaton decay and £. However, this does 




0.2 0.4 0.6 0.8 1 1.2 
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FIG. 3: One and two sigma contours for the 2D marginalized 
probability distribution, using CMB, 2dF and HST data, of 
the lepton degeneracy (£) and the curvaton density proportion 
at time of curvaton decay (r) for the case when lepton number 
is created by curvaton decay. 



not stop the BBN data strengthening the constraint on 
the isocurvature perturbation as the value of the curva- 
ton energy density at the time of curvaton decay is con- 
strained by the limits on the non-Gaussianity parameter 
from WMAP. So in this case, the neutrino isocurvature 
mode is less than about 5% (95% CI) of the adiabatic 
mode. Unless the BBN constraints on £ become tighter 
around zero, the case of the lepton number being cre- 
ated by curvaton decay would be unlikely if future ob- 
servations found / n i < —10, as this would imply a larger 
isocurvature perturbation than is observed. 
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